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Abstract 

Design methods for paraunitary matrices from complete orthogonal sets of idempotents and 
related matrix structures are presented. These include techniques for designing non-separable mul- 
tidimensional paraunitary matrices. Properties of the structures are obtained and proofs given. 
Paraunitary matrices play a central role in signal processing, in particular in the areas of filterbanks 
^SJ ' and wavelets. 

o 

1 Introduction 

A one-dimensional (ID) paraunitary matrix over C is a square matrix U(z) satisfying U(z)U* (z^ 1 ) = 1. 
Here * denotes complex conjugate transposed and 1 denotes the identity matrix of the size of U(z). In gen- 
eral a fc-dimensional (kD) paraunitary matrix over C is a matrix U (z) where z = (z\, Z2, ■ ■ ■ , Zk) is a vector 
of (commuting) variables {z\, z 2 , • • • , z^} such that U(z)U*(z~ 1 ) = 1 and z _1 = (z± , z% , . . . , z^ 1 ). 
Over fields other than C a paraunitary matrix is a matrix U(z) satisfying U(z)U(z~ 1 ) 1 - = 1. 
Paraunitary matrices are important in signal processing and in particular the concept of a paraunitary 
r/_ ■ matrix plays a fundamental role in the research area of multirate filterbanks and wavelets. In the 

polyphase domain, the synthesis matrix of an orthogonal filter bank is a paraunitary matrix; see for 
example [S]. 

Orthogonal filter banks may also be used to construct orthonormal wavelet bases [23 [T7] ; see also 
references in |14| . Paraunitary matrices over finite fields have been studied for their own interest and for 
, applications; see for example [16]. 

| Here general methods for constructing and designing such matrices from complete orthogonal sets of 

idempotents together with related matrix schemes are presented. This includes methods for designing 
I/*-) | non-separable multidimensional paraunitary matrices. Construction methods for complete orthogonal 

sets of idempotents are included. Group ring construction methods were the original motivation and 
from these more general methods evolved. A structure called the tangle of matrices is introduced; this 
may have independent interest. 

In certain cases specialising the variables of the paraunitary matrices allows the construction of 
series of regular real or complex Hadamard matrices. Walsh-Hadamard matrices, used extensively in the 
communications' areas, are examples of such regular Hadamard matrices. Complex Hadamard matrices 
arise in the study of operator algebras and in the theory of quantum computation. 

It is noted that the renowned building blocks for ID paraunitary matrices over C due to Belevitch and 
Vaidyanathan as described in [7] are constructed from W — {Fi,F 2 } where W is a complete orthogonal 
set of two idempotents in which F± has rank 1 and Fi has rank [n — 1) with n the size of the matrices 
under consideration. See section below for details on this. 

Connections between group rings, matrices and design of codes have been established in [5], [TU] and 
[1 lj : these are related but independent. 

Designing non-separable multidimensional paraunitary matrices is deemed difficult as there is no 
multidimensional factorisation theorem corresponding to the ID factorisation theorem of Belevitch and 
Vaidyanathan (}7;). For the finite impulse response (FIR) case there seems to be only a few examples 
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such as [IB] ■ See also [2] for background and further discussion. In [T3] a factorization of a subclass of 
2D paraunitary matrices is obtained; these though involve IIR (infinite impulse response) systems. 

In Section [5] results are obtained on the ranks of the idempotents and on the determinants of the 
paraunitary matrices formed. 

The concept of a pseudo-paraunitary matrix is introduced in Section [5] and construction methods for 
these are given. These may also be considered as FIR (finite impulse response) systems. 

2 Further Notation 

The book [3] is an excellent reference for background material on the algebraic structures used. 

Now F denotes a general field, R denotes a general ring, C denotes the complex numbers, K denotes 
the real numbers and Q denotes the rational numbers. Also ¥ q denotes the finite field of q elements, 
Rnxm denotes the set ofnxm matrices with coefficients from R and R[z] denote the polynomial ring 
with coefficients from R in commuting variables z = {zi, £2, • • • , z s ). Note that R[z] nxrn = R nxm [z]. 

Let R be a ring with identity li? = 1. (In general 1 will denote the identity of the system under 
consideration.) A complete family of orthogonal idempotents is a set {ei, e2, . . . , &k\ in R such that 

(i) 7^ and ef = e i; 1 < i < k; 

(ii) If i 7^ j then e^ey = 0; 

(iii) 1 = ei + e 2 + . . . + e k . 

The idempotent is said to be primitive if it cannot be written as = e i + e i where e i , e i are 
idempotents such that e i 7^ 0,e i 7^ and e i e i = 0. A set of idempotents is said to be primitive if each 
idempotent in the set is primitive. 

Various methods for constructing complete orthogonal sets of idempotents are derived below. Such 
sets always exist in FG, the group ring over a field F, when charF J(\G\. See [4] for properties of group 
rings and related definitions. These idempotent sets are related to the representation theory of FG. 

A mapping * : R — > R in which r 1— > r* , (r G R) is said to be an involution on R if and only if (i) 
r** = r, Vr G R, (ii) (a + 6)* = a* + b* , Va, b € R, and (iii) (06)* = b*a*, Vo, b € R. 

We shall be particularly interested in the case where * denotes complex conjugate transpose in the 
case of matrices over C and denotes transpose for matrices over other fields. Such a mapping * on group 
rings is also defined below. 

An element r € R is said to be symmetric (relative to *) if r* = r and a set of elements is said to be 
symmetric if each element in the set is symmetric. 

Q ® R denotes the tensor product of the matrices Q, R. 

As already noted A* is used to denote the complex conjugate transpose of a matrix A. Suppose R 
is a ring with involution *. Then * may be extended to matrices over R as follows. Let M G R n xm 
and define M* to be the matrix with each entry u of M replaced by u* . Then define M* — M* T . 
This matrix M* has size m x n. Let A(z) be a matrix with polynomial in variables z over some ring 
with involution *. Define A(z)* to be A*(z^ 1 ). When A is used for A(z) write A* to mean A(z)* . (In 
other words consider 'complex conjugate transposed' of a variable z to be this is consistent with 
group/group ring considerations.) 

Let R be a ring with involution *. For w(z) G R[z] define w(z)* = w*{z~ 1 ). Say w(z) is a paraunitary 
element in R[z] (relative to *) if and only if w(z)w* (z^ 1 ) = w(z)w(z)* = 1. 

Suppose K = (£>i, B2, ■ . ■ ,-Bfe) and L — {C\, C2, ■ ■ ■ , Ck) are rows of blocks of a matrix P where 
each block is of the same size. Then define the block inner product of K and L, written K ■ L, to be 
K ■ L = B\C* + B2C2 + ■ ■ ■ + BkCfr, This is to include the case when Bi, Cj are polynomial matrices 
and the C* are defined as above. 

3 Paraunitary elements 

The building methods using complete sets of orthogonal idempotents for the ID paraunitary matrices in 
this section are generalised later in section [5] below and following. The next section 2] considers methods 
for designing such complete sets of orthogonal idempotents. 
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Proposition 3.1 Let I = {ei, e 2 , . . . , e&} be a complete orthogonal set of idempotents in a ring R. 

k 

Define u(z) = ^""^ ±ejZ ti . Then u(z)u(z^ 1 ) = 1. 

i=l 

Proof: Since {ei, e 2 , . . . , e^} is a complete set of orthogonal idempotents, u(z)u(z^ 1 ) — e 2 +e 2 + . . -+e 2 = 
ei + e 2 + . . . + e k = 1. □ 

Corollary 3.1 J/ / is symmetric then u(z)u*(z~ 1 ) = 1. 

Thus u(z) is a paraunitary element when 7 is a symmetric orthogonal complete set of idempotents. 

It is not necessary to use primitive idempotents. Note also that if S = {e\, . . . , e^} is a complete set 
of orthogonal idempotents then {ej, ej}, i ^ j, may be replaced by {e^ + ej} in S and the result is (still) 
a complete set of orthogonal idempotents. This idea may be used to obtain real paraunitary matrices 
from (complex) complete orthogonal sets of idempotents in group rings. 

We single out the case R — F nxn for special mention. 

Proposition 3.2 Let {7i,7a, . . . , 7^} be a complete symmetric set of orthogonal idempotents in the ring 
F nxn of {n x n) matrices over F . Then W{z) = Y]j—i ±7^2** is a paraunitary ID n x n matrix over F 
where the U are non-negative integers. 

In the group ring case a paraunitary element in FG with \G\ — n gives a paraunitary matrix in F nxn 
via the embedding of FG into F nxn as given for example in 

Suppose {7i, 7 2 , ... , I k } is an orthogonal symmetric complete set of idempotents in F nxn and that P 
is a unitary matrix. Then also {P* IiP, P* I2P, ■ ■ ■ ,P*IkP} is a symmetric complete orthogonal set of 
idempotents in F nxn . 

For our purposes say a paraunitary matrix P is separable if it can be written in the form P = QR or 
P = Q ® R where Q, R are paraunitary with Q ^ 1,72 ^ 1; otherwise say P is non-separable. 

The following standard lemma is included for completeness and is not needed subsequently; the proof 
is omitted. 

Lemma 3.1 Suppose A(z) is a paraunitary matrix. Then A* (z) and A T (z) are paraunitary matrices. 
3.1 Modulus 1 

In Proposition 13. II the coefficients of the idempotents are ±1 times monomials. This can be extended in 
C to coefficients with modulus 1 times monomials. In R and fields of finite characteristic define a* = a 
and then ±1 are the only elements which satisfy aa* = a 2 = 1. 

Suppose {E\, E2- ■ ■ ■ ,Ek} is a complete symmetric orthogonal set of idempotents in F nxn . Define 
W{z) = a 1 E 1 z tl +a 2 E 2 z t2 + . . .+a k E k z tk and then W*^- 1 ) = a 1 *E 1 z- tl +a 2 *E 2 z- t2 +. . .+a k * E k z' tk . 
Here if a € C, then a* = a, the complex conjugate of a, and for other fields a* — a. Use \a\ 2 to mean 
00* for any field. 

Therefore W(z)W* (z^ 1 ) = W(z)W(z)* = \a 1 \ 2 E 1 + \a 2 \ 2 E 2 + ... + \a k \ 2 E k (**). 

Proposition 3.3 W(z) is a paraunitary matrix if and only if \a.i\ 2 = 1 for each i. 

Proof: If each \a t \ 2 = 1 then from (**) W(z)W*(z- 1 ) = 1. If on the other hand W(z)W*(z- r ) = 1 
then multiplying (**) through (on right) by Ei gives |ai| 2 i?i = Ei from which it follows that \cti\ 2 = 1. 
□ 

Thus Proposition 13. II may be generalised as follows: 

Proposition 3.4 Let {Ex, E 2 , . . . , E k } be a complete symmetric orthogonal set of idempotents and W{z) — 
ct\E\z l1 + a 2 E 2 z t2 + . . . + a k E k z tk , with tj > and \cej\ 2 — 1 for each j. Then W(z) is a paraunitary 
matrix. 

Now in C, \a\ 2 — 1 if and only if a = e ld for real 9 with i = \J — 1 and in R, |a| 2 = 1 if and only if 
a = ±1. In a field of characteristic p, \a\ 2 — a 2 = 1 if and only if a = 1 or a = —1 = p — 1. 
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As expected unitary matrices are built from complete symmetric orthogonal sets of matrices as per 
Proposition 13.41 

Proposition 3.5 U is a unitary nx n matrix over C if and only if U — a\V*Vi + a^ujjj v± + . . . + a n V*v n 
where {v\,V2, ■ ■ ■ ,v n } is an orthonormal basis for C„ and on G C, \on\ = 1, Vi. Further the on are the 
eigenvalues of U. 

Proof: Suppose U = aiV*V\ + azv^Vi + . . . + a n v* n v n with {v\, v^, . . . , v n } an orthonormal basis and 
\ati\ = 1. Then Uv* = onv* and so the on are the eigenvalues of U. It follows from Proposition 13.41 that 
U is unitary since {v*v\, v^V2, ■ ■ ■ , v^v n } is a complete symmetric orthogonal set of idempotents. 

Suppose then U is a unitary matrix. It is known that there exists a unitary matrix P such that 

U = P*DP where D is diagonal with entries of modulus 1. Then P = I I where {v\,V2, ■ ■ ■ ,v n } is 

an orthonormal basis (of row vectors) for C n and D = diag(ai, Q!2, . . . , a n ) with \on\ = 1 and the on are 
the eigenvalues of U. Then 

U = P*DP 

(Ql o ... o 
Q 2 ... 
: : : : 
6 6 dr. 

= (aivl,a2V^...,a n v^) 
= aiv*Vi + a2V%V2 + ■ ■ ■ + a n v*v n . 

□ 

Thus unitary matrices are generated by complete symmetric orthogonal sets of idempotents formed 
from the diagonalising unitary matrix. Notice that the on are the eigenvalues of U. 

For example consider the real orthogonal/unitary matrix U = ( _° shfe co" e ) ■ This has eigenvalues 
e i , e~™ and P = A= y?- ^ ! ) is a diagonalising unitary matrix. Take the rows «i = -^^( — 1, — z), V2 — 
~7|(*)1) °f P an d consider the complete orthogonal symmetric set of idempotents {Pi = V*v\ — 

i(i^).ft=^ = i(i.i)}- 

Then applying Propositon 1331 gives U = e ie P x + e~ ie P 2 = \e ie (] ~*) + |e" ie (ij 1), which may be 
checked independently. 

3.2 Products 

A product of paraunitary matrices and the tensor product of paraunitary matrices are also paraunitary 
matrices. Thus further paraunitary matrices may be designed using these products from those already 
constructed . 





4 Complete orthogonal sets of idempotents 

Paraunitary matrices are designed from complete symmetric sets of orthogonal idempotents in section [3] 
and also in later sections. Here we concentrate on how such sets may be constructed. 

4.1 Systems from orthonormal bases 

Let V — F n . Assume F n has an inner product so that the notion of orthonormal basis exists in V 
and its subspaces. In R" and C™ the inner product is vu* for row vectors v,u where * denotes complex 
conjugate transpose; in K", w* = w T , the transpose of w. 
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Suppose now V = V\ ffi V2 ffi • ■ • ffi 14 is any direct decomposition of V. Let Pi denote the projection 
of V to Vi. Then Pj is a linear transformation on V and (i) 1 = Pi + P2 + . . . + P&; (ii) Pf — Pi] (hi) 
/'/'., 0. <•'./. 

Thus {Pi, P 2 , . . . , Pfc} is complete orthogonal set of idempotents. If each Pj is an orthogonal projection 
then this set is a complete symmetric orthogonal set of idempotents. 

The matrix of Pi may be obtained as follows when Pj is an orthogonal projection. Let {w\ , W2 , • • • , w s } 
be an orthonormal basis for Vi and consider w £V. Then w = Vi + w where w £ V\ © V2 ffi ■ ■ ■ Vi . . . © 14 
and Vi £ Vi. Here Vi means omitting that term. Then Pj : V Vi is given by hi Uj. Now 

= a\Wi + ot2W2 + . . . + ctkWk- Take the inner product with Wj to get ctj = ViW* = ww*. Hence 
Pi : w >->■ w(w{wi + u> 2 u> 2 + . . . + w*w s ). Thus the matrix of Pj is rofwi + W2W2 + • • • + w*w s . 

On the other hand suppose {Pi, P2, . . . , P&} is a complete symmetric orthogonal set of idempotents 
in F n xn- Then p defines a linear map V — > V by Pj : u uPj. Let Vi denote the image of Pj. Then it 
is easy to check that V = Vi © V 2 © . • • ffi 14- 

The case when each Vi has dimension 1 is worth looking at separately. Suppose {01, o 2 , . . . , o„} is 
an orthonormal basis for F n . Such bases come up naturally in unitary matrices. Let Pj denote the 
projection of P™ to the space generated by o^. Then P = {Pi, P2, . . . , P n } is an orthogonal symmetric 
complete set of idempotents in the space of linear transformations of P" . It is easy to obtain the matrices 
of Pi . The matrices Pj may be combined and the resulting set is (still) a complete symmetric orthogonal 
sets of idempotents. For example (Pj + Pj) (i ^ j) is still idempotent and is the projection of F n to the 
space generated by {oj,Oj}; replace {Pj,Pj} by (Pj + Pj) in P and the new set is (still) an orthogonal 
symmetric complete set of idempotents. Then rank(Pj + Pj) = rank(Pj) + rank(P 3 ) also - see Lemma 
ETT1 below. 

For example {v% = |(2, 1,2),U2 == ^ (1, 2, — 2),V3 — |(2, — 2, — 1)} is an orthonormal basis for R 3 . 
The projection matrices are respectively Pi = Vi Vi = ^ ( 2 1 2 J , P 2 = V2 v 2 = | ( 2 4 -4 J , P 3 = 

«3 T ^ = ^( = 4 4 4 4 f). 

Thus {Pi,p2,P3} is a complete symmetric orthogonal set of idempotents and each Pj has rank 1. 
Set P 2 = P 2 + P3 and then {Pi,p2} is a complete symmetric orthogonal set of idempotents also and 
rank(P 2 ) = 2. 

Note that the inner product in C™ is vu* for row vectors v,u where * denotes complex conjugate 
transposed. For example h !)> 7^ (i, 1)} is an orthonormal basis for C 2 . Projecting then gives the 

complete orthogonal symmetric set of idempotents {Pi = | ( \ ~£ ) , P 2 = \ ( 2j \ )}• 

4.2 Orthogonal idempotents systems from unitary /paraunitary 

Let U be a unitary or paraunitary nx n matrix in variables z say over R. Then the rows {vi, V2, ■ ■ • , w n } 
of J7 satisfy WjW* = 1, UjW* = 0, i 7^ j. 

Define Pj = v,*Uj for i — 1, 2, . . . , n. Then Pj is an n x n matrix of rank 1. 

Proposition 4.1 {Pi, P2, • • • , P n } * s a complete symmetric orthogonal set of idempotents in P„ x „[z,z -1 ]. 

Proof: It is easy to check that P* = Pj, PjPj = Pj, PjP, = 0, i 7^ j. It is necessary to show that the set 
is complete. 

Let A = Pi + P 2 + . . . + P n . Note that P t v* = v*,P z v* = 0, i ^ j. Then Av* = v*. Thus A has n 
linearly independent eigenvectors corresponding to the eigenvalue 1. Hence A — I n . □ 

4.2.1 Diagonals 

In R n xn let Pjj denote the matrix with 1 = 1^ on the (diagonal) (i,i) position and elsewhere. Then 
W = {Pn, P22, • • • j Enn} is a complete symmetric orthogonal set of idempotent matrices. This is a 
special case of section 14.11 but is worth mentioning separately; paraunitary matrices have been designed 
from W which, although not generally useful in themselves directly, may be combined with other designed 
paraunitary matrices with which they do not commute in general. 
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4.3 Group rings 



Group rings are a neat way with which to obtain complete orthogonal symmetric sets of idempotents. 
These systems have nice structures from which properties of the paraunitary matrices designed may be 
deduced. Let w = X^eG a g9 ^ e an element in the group ring FG and W denotes the matrix of w as 
defined in [IT] and [9]. This matrix W depends on the listing of the elements of G and relative to this 
listing tj> : w i-> W is an embedding of FG into the ring ofnxn matrices, F nxn , over F where n = \G\. 
The transpose, w T , of w is w T = X s gg oigg^ 1 . Note that the matrix of w T is then W T . 

Over C define w* — ^2 ge G^99 where bar denotes complex conjugate. Note then that for a group 
ring element w with corresponding matrix W the matrix of w* is indeed W* . 

Say an element w(z) G FG[z] is a paraunitary group ring element if and only if w(z)w* (z _1 ) = 1 and 
this happens if and only if the corresponding W{z) G F nxn [z] is a paraunitary matrix (where n = \G\). 
The W(z) obtained from w(z) in this case is termed a group ring paraunitary matrix. 

Group rings are a rich source of complete sets of orthogonal idempotents and group rings have a rich 
structure within which properties of the paraunitary matrices so designed may be obtained. 

The theory brings representation theory and character theory in group rings into play. The orthogonal 
idempotents are obtained from the conjugacy classes and character tables, see e.g. [2|. The orthogonal 
sets of idempotents depend on the field under consideration and classes of paraunitary matrices over 
different fields such as Q, R or finite fields are also obtainable. 

The primitive central idempotents of the complex group algebra CG are given by e(x) = X g eG x(9 1 )9 
where \ runs through the irreducible (complex) characters \ °f see [4], Theorem 5.1.11 page 185, 
where the are expressed as = ^tqt X s <eg Xi(d~ 1 )9- 

The idempotents from group rings are automatically symmetric. 

Theorem 4.1 For the group idempotents e-i, e* = e\. 

Proof: This is a matter of showing that the coefficients g and g _1 in each ei are complex conjugates of 
one another. But this is immediate as it is well-known that xG? _1 ) — x(fiO) an d thus the result follows 
from the expression for given above. □ 
Let Ei denote the matrix of as per an embedding of the group ring into the ring of matrices as for 
example in |11) . 

Corollary 4.1 Let {ei, e^, ■ ■ ■ , e^} be a complete set of orthogonal idempotents in a group ring and define 
fe 

U(z) — zLEiZ** where the U are non-negative integers. Then U(z) is a paraunitary matrix. 

i=l 

Corollary 4.2 Let {ei, e2, . . . , e^} be a complete set of orthogonal idempotents in a group ring over C 
fc 

and define U(z) = ^^aiEiZ li where the U are non-negative integers and \<Xi\ = 1. Then U{z) is a 
i=i 

paraunitary matrix. 

The formula for the {ei} as given above (taken from [3]) may be used to construct complete orthog- 
onal sets of idempotents. The Computer Algebra packages GAP and Magma can construct character 
tables and conjugacy classes from which complete sets of orthogonal idempotents in group rings may be 
obtained. The literature contains other numerous methods for finding complete (symmetric) orthogonal 
sets of idempotents in group rings. 

In general the paraunitary matrices designed using orthogonal sets of idempotents in the group ring 
over C have complex coefficients but specialising and combining idempotents allows the design so that 
the coefficients may be in R, the real numbers, or in Q, the rational numbers. When the group ring of 
the symmetric group S n is used the paraunitary matrices derived by these methods all have coefficients 
automatically in Q and when the group ring of dihedral group Z?2n is used the coefficients are in R. In 
general idempotents occur in complex conjugate pairs and these may be combined to give real coefficients 
resulting in paraunitary matrices with real coefficients. 

Most of the results hold in the case when the characteristic of F does not divide the order of G; in 
this case this means that the characteristic of F does not divide the size n of the (n x n) matrices under 
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consideration. In these cases also it may be necessary to extend the field to include roots of certain 
polynomials. 



4.4 Tensor products 

It is easy to check that the tensor product of paraunitary matrices is also a paraunitary matrix. If 
P = QR, S = TV then P®S = (QT) ® (RS) when the products QT, RS can be formed. 

Complete orthogonal sets of idempotents may be designed using products of these sets. Suppose 
{eo, e2, . . . , efe} is a complete orthogonal set of idempotents in F nxn and {fo, fi, ■ ■ ■ , f s } is a complete 
orthogonal set of idempotents in Ff. X k- Then {e* ® fj | < i < k, 1 < j < s} is a complete orthogonal set 
of idempotents in F n k Xn k- Here ® denotes tensor product. If both {eo, e2, • • • , efc} and {/o, /i, . . . , / s } 
are symmetric then so is the resulting tensor product set. The details are omitted. 

If {ei | 1 < i < k} and {fj | 1 < j < s} are complete orthogonal sets of idempotents within group 
rings FG, FH respectively then {eifj |l<z</c,l<j<s}isa complete orthogonal set of matrices in 
F(G x H). Suppose e^ i-> Ei, fj i-> Fj gives an embedding into matrices, then eifj H> Ei® Fj gives an 
embedding into F{G x H); this may be deduced from [TT] and details are omitted. 



4.5 Examples of paraunitary matrix from orthonormal bases 

i/424\ 1 / 1 2 — 2\ 

1. The complete orthogonal symmetric systems of idempotents Pi = | (212) , P2 — ^ ( 2 4 -4 J , _P 3 

. / 4 -4 -2 \ 

j 1-4 4 2 1 were obtained in section l4.1l Then W(z) = Piz +P2Z+P3Z is a paraunitary matrix. 

2. Let z = e l9 in in 1. gives a unitary matrix, T say. The rows of T form an orthonormal basis 
for C 3 . These rows may then be used to form a complete symmetric orthogonal set of idempotents 
from which paraunitary matrices may be constructed. This process could be continued. 

3. In CC3, where C3 is the cyclic group of order 3, the orthogonal complete set of idempotents formed 

are G\ — h f 1 1 1 ^ ,02 = I I ^ m 1 u | , Qs = h ( u l^ 2 ) where uj is a primitive 3rd root of 
3 Vi 1 r) 3 U" 2 1/ V" 2 u 1 / 

unity. 

Then Q(z) = Q% + Q2Z 3 + Q3Z 2 is a paraunitary matrix. 

4. Give values of modulus 1 to z in Q(z) above and get a unitary matrix R. Use the rows of R to 
form a further complete symmetric set of idempotents from which paraunitary matrices may be 
formed. 

5. Combine Q(z) in 3. with W(z) in 1. to give for example the paraunitary matrix Q(z)W(z)Q(z). 

We give some examples from orthogonal sets of idempotents derived from group rings. The group ring idea 
is used later as a prototype in which to extend the method for the design of non-separable paraunitary 
matrices. The complete orthogonal sets of idempotents obtained from group rings are automatically 
symmetric as noted in Theorem 14. II 

Recall that circ(ao, 01, . . . , a n —i) denotes the circulant n x n matrix with first row (ao, fflj., . . . , a n -i)- 

Consider CC„ where C n is a cyclic group of order n. 

1. When n — 2 the (primitive) orthogonal set of idempotents consists of {eo = l/2(l + g), e\ = 1/2(1 — 
<?)}, where g generates C%. Thus paraunitary matrices may be formed from Eq = | f j | J and 

Ex = I ^ ^ ^ J giving for example - (Looks familiar?) 

2. These may be combined with paraunitary matrices formed from En = \ \ ^ ) ,^22 ' 



Note that En, E22 do not commute with Eq,E\. For example the following is a paraunitary matrix: 
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fl 0\lfz + z 2 z-z 2 \(z 2 0\l/z 2 + z 3 z 2 -z 3 \ 
[O z) 2 \z - z 2 z + z 2 ) \ z 3 ) 2 \z 2 - z 3 z 2 + z 3 ) 

They may also be combined with paraunitary matrices formed from orthonormal bases as in section 
Osuch as {P x = ,P 2 = | (2 i *)}. 

The determinant of these matrices which are powers of z may be obtained from Theorem |H2] below. 

3. The primitive orthogonal idempotents for a cyclic group are related to the Fourier Matrix. 

4. In CC 4 , for example, the orthogonal primitive idempotents are e 4 = 4(1 + a + a 2 + a 3 ),e2 = 
i(l + w 2 a 2 + w 3 a 3 ), e 3 = - a + a 2 - a 3 ), e 4 = j(l +uj 3 a + u 2 a 2 + ujo 3 ) from which 4x4 
paraunitary matrices may be constructed. Here w is a primitive 4 th root of unity and in this case 
uj 2 = -1. " 

Notice that e, = e^ 1 as could be deduced from Theorem 14. II 

5. Combine the e^ to get real sets of orthogonal idempotents. Note that it is simply enough to combine 
the conjugacy classes of g and g" 1 . In this case then we get 

ei = e\ = 4(1 + a + a 2 + a 3 ), e 2 = e 2 + e 4 = i(l — a 2 ), = e 3 — 4(1 — a + a 2 — a 3 ), which can 
then be used to construct real paraunitary 4x4 matrices. 

6. Using C2 x C2 gives different paraunitary matrices. Here the set of primitive orthogonal idempotents 
consists of fi = j(l + a + b + ab), / 2 = — a + b — ab), = ^(1 — a — b + ab), / 4 = |(1 + a — 6 — a6) 
and the paraunitary matrices derived are all real. 

7. The paraunitary matrices produced from C4 from C 2 x C 2 and from En, E22, E33, E44 may then 
be combined to produce further (4 x 4) paraunitary matrices. So for example the following 4x4 
is a paraunitary matrix: 

{Ex + E 2 z + E 3 z 3 + E 4 z 2 ){En + E 22 z + E 33 z 3 + E 4i z 2 )(Fiz + F 2 z 2 + F 3 z 3 + F iZ 2 ) 

Again the determinant of the matrix may be obtained from Theorem 19.21 

The Ei,Fj are derived from the fj (as per [TT]) so for example E2 = \ ciic(l,u,u 2 ,uj 3 ), F 3 = 
/ 1 -x -1 1 \ 

X / -X X X -X ) 
4 1 -X X X -X I ■ 

\ X -1-1 1 / 

4.6 Get real 

By combining complex conjugate idempotents in a complete orthogonal sets of complex idempotents, 
real paraunitary matrices may be obtained. We illustrate this with an example. 

Suppose {eo, e±, e 2 , e 3 , e 4 , es} is the complete set of primitive idempotents in CCq. Here then — 
|(1 + u l g + uj 2l g 2 + Lo 3l g 3 + uj^g 4 + D 5l g 5 ) where uj — e 2l7r / 6 is a primitive 6 th root of unity and Cq is 
generated by g. 

Then eJJ = e ,eT = e 5 ,ei = e 4 ,e^ = e 3 . 

Let = 27r/6. Note that co&{6) — cos(50), cos(26') = cos(4(9). Now combine ei with es and e 2 with 
e 4 to get ei = |(1 + cos(0)g + cos(26)g 2 + cos(3% 3 + cos(46»)g 4 + cos(5% 5 ) and e 2 = |(1 +cos(20)g + 
cos(29)g 2 ) + cos(28)g 3 + cos(26)g 4 + cos(28)g 5 . This gives the real orthogonal complete set of idempotents 
{eo, ei, e2, 63} from which real paraunitary matrices may be constructed. The ranks of the idempotents 
and determinants of the paraunitary matrices formed may be deduced from Lemma 19.11 and Theorem 



4.7 Symmetric, dihedral groups 

Let £) 2n denote the dihedral group of order 2n. As every element in _D 2n is conjugate to its inverse, the 
complex characters of -D 2ra are real. Thus the paraunitary matrices obtained directly from the complete 
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orthogonal set of idcmpotents in C£>2n have real coefficients. The characters D 2n are contained in an 
extension of Q of degree <f>(n)/2 and this is Q only for 2n < 6. 

Let S n denote the symmetric group of order n. Representations and orthogonal idempotents of 
the symmetric group are known; see for example [5]. The characters of S n are rational and thus the 
paraunitary matrices produced directly from the complete orthogonal set of idempotents in CS n have 
rational coefficients. 

The paraunitary matrices formed from different group rings (with same size group) may be combined 
to form further paraunitary matrices; these in general will not commute. 

We present an example here from S3, the symmetric group on 3 letters. (Note that S3 = Dg.) 
Now 53 = {1, (1, 2), (1, 3), (2, 3), (1, 2, 3), (1, 3, 2)} where these are cycles. We also use this listing of 
S3 when constructing matrices. 

There are three conjugacy classes: K 1 = {1}; K 2 = {(1, 2), (1, 3)}, (2, 3); K 3 = {(1, 2, 3), (1, 3, 2)}. 
Define 

ei = 1 + (1,2) + (1, 3) + (2, 3) + (1, 2, 3) + (1, 3, 2), 
e 2 = 1- {(1,2) + (1,3) + (2,3)} + (1,2,3) + (1,3,2), 
e 3 = 2- {(1,2,3) + (1,3,2)}, 

and ei = \&i;e 2 = \e2\e3 = ^£3. Then {ei, 62,63} form a complete orthogonal set of idempotents 
and may be used to construct paraunitary matrices. 

/ 1 (12) (13) (23) (123) (132) 



The G- matrix of S3 (see [IT]) is 



(12) 1 (132) (123) (23) (13) 

(13) (123) 1 (132) (12) (23) 
(23) (132) (123) 1 (13) (12) 
(132) (23) (12) (13) 1 (123) 

\(123) (13) (23) (21) (132) 1 

Thus the matrices of e%, e 2 , e3 are respectively 




1-1-1-11 i 
-li i l-i-i 
-li i l-i-i 
-li i l-i-i 
1-1-1-11 i 
l-i-i-ii i 



, E3 



2 











-1 


-1 





2 


-1 


-1 











-1 


2 


-1 











-1 


-1 


2 








-1 











2 


-1 


-1 











-1 


2 



Note that E\,E2 have rankl and that E3 has rank 4; the proof for the ranks of these Ei in general 
is contained in Lemma |9. II 

Thus for example the following are paraunitary matrices: 
E 2 + Eiz + E 3 z 2 , E3 + {Ex + E 2 )z, Ex+E 2 + E 3 z 2 . 

The paraunitary matrices formed from these idempotent matrices may then be combined with pa- 
raunitary matrices formed using complete orthogonal idempotents obtained from CCq and ones using 
{En, E22, E33, E44, E 55 , Eqq\. 



Let {/1, . . . , fe} be the orthogonal idempotents from CCq. Let wi = J2i=i E{Z li , w 2 = J2i=i 



Si=i EiZ Li . (Fi is the matrix of fi ) Then products of w\, W2,W3 are paraunitary matrices. Note that 
the Wi do not commute. 



4.8 Finite fields 



Here we consider constructing examples of complete symmetric sets of orthogonal idempotents over finite 
fields. 

Suppose {vi, «2, . . . , Vk} is a orthogonal basis for F k under (it, v) = uv T . Thus ViVj T = for i + 1 j. 
Suppose also (vi,Vi) = ti + 1 0. Define Pi = t~ 1 Vi T Vi which is a k x k matrix. Then PiPi T = 
t~ 1 v i T v i t~ 1 v l T v i = t~ 1 v T v i = Pi and PiPj T = t~ 1 v i T v i tj 1 v j T v : j = for i + 1 j. 

k 

It also follows that Pj =1. To see this consider A = Pi + P 2 + . . . + Pk- Then v, L A = UjPj = vi as 
i=i 

hi\ 

1>2 

ViPj = for i 7^ j and ViPi — Vit~ 1 v{ T Vi = Vi. Hence ViA — Vi. Let Q = 



Then Q is non-singular 
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as {vi, V2, ■ ■ ■ , Vk} is linearly independent. Then QA — Q and hence A = 1. 
Note that in the above we do not need to take the square root of elements. 

Another way could be to construct such sets over Q and when the denominators do not involve a 
prime dividing the order of the field it is then possible to derive complete symmetric orthogonal sets of 
idcmpotents over the finite field. 

For example the complete orthogonal symmetric systems of idempotents Pi = § ( 2 l 2 J , P 2 = 



2 -2\ . / 4 -4 -2 \ 

4 -4 ) , P 3 = i ( -4 4 2 ) were obtained in section I4TU 



-2-44/ J V -2 2 1 

Over a field of characteristic 2 these come to the trivial set {(010), (000), (000)} of symmetric 

L \o 0/ Vo 0/ Vo 1/ J 

complete orthogonal set of idempotents. 

Over the field F5 of 5 elements they become (note that here 9~ 1 =4): {f^|^))f|||)'f|||)}- 

This is a complete symmetric orthogonal set of idempotents in F5 which may be checked indepen- 
dently. 

The following are complete symmetric orthogonal sets of idempotents over F7: 

r /212\ /416\ /2 5 6\, r /6 5 6\ / 5 2 5 \ / 4 3 \ 
{141,125 , 521}, { 535 , 252 , 000 }. 
L V212/'V652/'V614/ J L V656/ V525/ V304/ J 

These different sets may be used to construct paraunitary matrices over F7 and in a later section are 
used to show how to construct as an example a non-separable paraunitary matrix over a finite field. 



4.9 ID building blocks 



The great factorisation theorem of Belevitch and Vaidyanathan, see [7] (pp. 302-322), is that matrices 
of the form H(z) = 1 — vv* + zvv*, where v is any unit column vector (v*v = 1), are the building blocks 
for ID paraunitary matrices over C. 

Consider F± = vv* where v is a unit column vector and so v*v — 1. Thus F±F± — vv*vv* — vv* = F\ 
and so F± is an idempotent. Hence {F± — vv* , F2 — 1 — F\ — 1 — vv*} is a complete symmetric orthogonal 
set of these (two) idempotents with rank.Fi = 1 and ranki<2 = (n — 1) where the matrices have size 
n x n; see Theorem 19.11 below for rank result. Then H(z) = F2 + zF\ and hence the paraunitary ID 
matrices are built from complete symmetric orthogonal sets of two idempotents, one of which has rank 1 
and the other has rank(n — 1). 



Proposition 4.2 Let F be a field in which every element has a square root. Suppose also an involution 
* is defined on the set of matrices over F. Then P is a symmetric (with respect to * ) idempotent 0/ 
rankl in F nxn if and only if P — vv* where v is a column vector such that v*v = 1. 



(Note that 'symmetric with respect to *' in the case of C in which * denotes complex conjugate 
transposed is normally termed 'Hcrmitian'.) 

Proof: If P = vv* with v*v = 1 then P is a symmetric idempotent of rank 1. 

Suppose P is a symmetric idempotent of rankl in F nxn . Since P has rank 1 each row is a multiple 
of any non-zero row. Suppose the first row is non-zero and that the first entry of this row is non-zero. 
Proofs for other cases are similar. Since P is symmetric it has the form 

/h b 2 ... b n \ 
I b% b 2 b*/b 1 ... b n b* 2 /b l 1 
63 b 2 b^/bi ... b n b^/bi 



P = 



\K b 2 bi/b t b„b*/bi / 

with b\ =b\. 

Since P is idempotent it follows that b\ + \b 2 \ 2 + ■ . ■ + \b n \ 2 =b\. 

Let v = ^(bi,b 2 , ... , b n )*. Then v*v = + |6 2 | 2 + . . . + \b n \ 2 ) = 1 



and vv* = -r- 



/ bihi 6162 ■■■ bib n 
b^bi b^b 2 ... b 2 b n 



\K b i b i b 2 b ib„ 



P. 



□ 
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It is necessary that square roots exist in the field. For example P = ( 2 \ ) over F3 is a symmetric 
idempotent matrix of rank 1 but cannot be written in the form vv 1 -; however 2 does not have a square 
root in F3. Note that P\ = 1 — -P = ( 2 2) an d that {P, Pi} is a complete orthogonal set of idempotents 
in (F 3 ) 2X 2- 

Over F3 the following complete symmetric sets of idempotents are the building blocks for 2 x 2 
matrices: {(12)5(22)}) {(00)5(01)}- Thus the paraunitary matrices 2x2 matrices over F3 are built 
from these sets using Proposition 13.41 and products. These sets are not of Belevitch and Vaidyanathan 
form and their result does not apply here. 

The ID building block result of Belevitch and Vaidyanathan cannot be extended to multidimensions. 

Group ring ID paraunitary matrices and other ID paraunitary matrices over C constructed here can 
in theory then be obtained from this characterisation. Group rings have special features and paraunitary 
matrices from these have nice structures. 

5 Multidimensional 

5.1 kD with idempotents 

Proposition 5.1 Let {E\, E2, ■ ■ ■ , E t } be a complete symmetric orthogonal set of idempotents and de- 
fine products of non-negative powers of variables by Wi(z) = ± Ilj=i z j 3 f or * = 1,2, where 

z = (zi, Z2, . . . , Zk) and tij are non-negative integers. Define W(z) = Y^i=i w i( z )Ei- Then W(z) is 
a k-dimensional paraunitary element. 

Proof: Since {Ex,E 2 , ■ ■ ■ ,E t } is a complete symmetric orthogonal set of idempotents, W(z)W(z)* = 
E 2 +E 2 + ...+E 2 =E 1 +E 2 + ... + E t = l. □ 
As before for a € C define a* = a, the complex conjugate of a, and for other fields define a* = a. 
Define \a\ 2 — aa* in all cases. 

Proposition 5.2 Let {E±, E2, ■ ■ ■ , E t } be a complete symmetric orthogonal set of idempotents and define 
products of non-negative powers of the variables by Wi(z) — cti Ylj—% Zj' 3 f or i = 1,2, ... ,t where z = 

(z\, 22, ... 1 Zk) and tij are non-negative integers and \cti\ 2 = 1. Define W(z) = J2l=i w i( z )Ei and then 
W(z) is a paraunitary matrix. 

The W(z) so formed can be combined using products of matrices, or tensor products of matrices 
when appropriate, to form further paraunitary matrices. 

Such paraunitary matrices formed from Propositions l5.1 1and l5.2l can however be shown to be separable 
but have uses of their own and will be used later to form (constituents of) non-separable paraunitary 
matrices. 

5.2 Examples of 2D paraunitary 

Recall that if u is a group ring element of FG with \G\ — n then U (capital letter equivalent) denotes 
the matrix of u under the embedding of FG into the ring F nXn of n X n matrices over F, see [TT] . 

Let {eo, ei, e 2 } be a (primitive) orthogonal complete set of idempotents in CC3, and thus {Eq, E\, E 2 } 
is an orthogonal complete symmetric set of idempotents in C3 X 3- Let {fo, fi, f 2 , fs\ be an orthogonal 
complete set of idempotents in CC4. Define u(z, y) = (eo + e\z + e 2 z 2 )fo + (eo + e%z + e 2 z 2 )f\y + (eo + 
e\z + e 2 z 2 )f 2 y 2 + (eo + e\z + e 2 z 2 )f^y i and let U(z, y) be obtained from u(z, y) by replacing each e^ by 
Ei and each by Fj. 

Then U{z,y) is a paraunitary matrix. However here u{z,y) = (eo + eiz + e 2 z 2 )(fo + fiy + f 2 y 2 + fsy 3 ) 
and so U(z, y) is separable. 

Let u(z,y) = (e Q + e 1 z + e 2 z 2 )f + (e + eiz 2 + e 2 z)fiy+(eoz + ei+e 2 z 2 )f 2 y 2 + (e z + eiz 3 + e 2 z 2 )f 3 y 3 . 
Then U(z, y) is a paraunitary matrix. 
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6 Matrices of idempotents 

Paraunitary matrices may also be constructed from matrices with blocks of complete orthogonal sets of 
idempotents. 

Consider the following example. Let Eq — ^ ( } } ) , E\ — | ( ~£ ) . 



( xE yE ± \ _ i 
I zEi tE J 2 





fx 


X 


y 


-y\ 


1 


X 


X 


-y 


y 


2 


z 


—z 


t 


t 






z 


t 


t J 



Define W 

Then WW* = I4 as {Eq,Ei} is an orthogonal symmetric complete set of idempotents. However W 
is separable as 

fxEo+Ex \{Eo yE 1 \ 

\ tEo + Ek) \zE 1 E J- [ > 
and each of the matrices on the right in (**) is separable into ID paraunitary matrices. 
Here if we let x = 1 = t then (**) is a trivial product as the first matrix on the right is then the 

identity. If further y = 1 = t this produces the matrix H = ^ I \ _} 1 ^ \ I which is a common matrix 

V -1 1 1 1 J 

used in quantum theory as non-separable. (This matrix H with fraction omitted is a Hadamard regular 
matrix.) Thus non-separability of a paraunitary matrix is a stronger condition by comparison. 

Let q = Hi D.Q^Kir)- ^=(aa)- 

Then W is a paraunitary matrix. Now letting the variables have complex values of modulus 1 gives 

/ 1 i 1—2 

rise to complex Hadamard regular matrices as for example [ ~i ^ \ \ 

V i 1 -i 1 

Let {Eq, Ei, E2} be an orthogonal symmetric complete set of idempotents in F3 X 3. 

(xE a yE x zE^ 
Define W = pE 2 qE a rE\ 

\sE x tE 2 vE 0/ 

The variables of W are x, y, z,p, q, r, s, t, v which need not necessarily be distinct. 
Then WW* = I g . 

For example in section T4. II the following complete set of symmetric idempotents was constructed in 

Q 3 x 3 : 

rp -1 /4 2 4\ T / 1 2 -2\ _ 1 / 4 -4-2 

Pi = Vl T Vl = i (2 1 2j ,P 2 = V 2 T V 2 = § (_2 2 _^ -4 j ,P 3 = v^v 3 = i (-4 4 2 

/xPi yP 2 zP 3 \ 
Then W = pP% qP\ rP 2 is a paraunitary matrix. 
\ S P 2 tP 3 vPj 



6.1 General construction 

Let {Eq, Ei, . . . , Ek} be a complete symmetric orthogonal set of idempotents in F nxn . Arrange these 
into a k x k block matrix of with each row of blocks containing one of the blocks {Eq, E\, . . . , Ek} exactly 
once. Now attach monomials to each Ef, the same monomial need not be used with each Ei that appears. 
Let W be the resulting matrix. 

Theorem 6.1 W is a paraunitary matrix. 

Proof: Take the block inner product of two different rows of blocks. The Ei are orthogonal to one 
another so the result is 0. Take the block inner product of any row of blocks with itself. This gives 
E\ + El + . . . + El = Ei + E 2 + . . . + E k = 1(= /„). Hence WW* = 1(= I nk ). □ 
The W is a paraunitary matrix in the union of the variables of the monomials. 

The condition that each row and column block contains each Ei once and once only can be obtained 
by using the group ring matrix of any group of order k; see for example So for example the Ei could 
be arranged as a circulant block of matrices. Different arrangements will in general give inequivalent 
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paraunitary matrices. 

(Modifications in the construction of W by attaching elements of modulus 1 as coefficients will give 
paraunitary matrices but further conditions are necessary on these elements. These modifications are not 
included here.) 

These are nice constructions for paraunitary matrices but can be shown to be separable. However they 
have uses in themselves and will prove useful later as parts of constructions of non-separable paraunitary 
matrices. They may also be used to construct special and regular types of Hadamard real and complex 
matrices. These are illustrated in the following sections 16.21 and [6.31 



6.2 Monomials and Hadamard regular matrices 

Although the matrices of idempotents as constructed in sections IBI and ROl produce separable paraunitary 
matrices these can be useful structures in themselves; they may also be used in certain cases to produce 
'regular' Hadamard matrices. (Walsh-Hadamard matrices are regular Hadamard matrices which have 
been used extensively in communications' theory.) If in a paraunitary matrix the entries are ±1 times 
monomials in the variables then substituting ±1 for the variable gives a Hadamard matrix. If in a 
paraunitary matrix the entries are to times monomials where w is a complex number of modulus 1 then 
substituting each variable by a complex number of modulus 1 gives a complex Hadamard matrix. 

For example use Pq = circ(l, 1, l),Pi = circ(l,w,a; 2 ),P2 = circ(l, uj 2 , w) where ui is a primitive third 

/xP yPx zP 2 \ 
root of unity gives the following matrix: zP 2 xPq yP\ 

\yP x zP 2 xP a ) 

Substituting a complex number of modulus 1 for each of x, y, z gives a complex Hadamard matrix. 
Butson-type Hadamard matrices H(q,n) are complex Hadamard n x n matrices with entries which 
are q th roots of unity. 

With the above example substituting a third root of unity for the variables gives a Hadamard H(3, 9) 
matrix, that is a a matrix H with entries which are third roots of unity so that HH* — 9. These matrices 
could then be used to produce Hadamard 7J(3,36) matrices. 

This can be extended to q x q matrices using the complete orthogonal set of idempotents for the cyclic 
group of order q; this will involve q th roots of unity and is related to the representation theory of the 
finite cyclic group. From this Hadamard H(q, q 2 ) matrices can be produced and from these Hadamard 
H(q, (2q) 2 ) matrices can be produced and so on. 



6.3 Mixing 



It is noted that interchanging rows and/or columns in a paraunitary matrix results in a paraunitary 
matrix and thus interchanging blocks of rows and/or columns results in a paraunitary matrix. When 
using complete orthogonal sets of idempotents the blocks of row and/or columns of idempotents can be 
interchanged in the construction stage. The resulting paraunitary matrices can take a regular form. In 
certain cases the variables can be specialised to form regular Hadamard matrices. Here examples are 
given but details of the constructions are omitted. 

Let {Eq, Ei} be a complete symmetric orthogonal set of idempotents in C nxn . For example in C 2x2 
these could be E = \ ( \ \ ) , E x = § ( _\ ) . 



Define 



W = \ 



( 



V 



xE 

-uEx 

zE 

-wEi 

xE 

uEx 

-zE 

-wE\ 



vE 
-tE 1 
pE 

yEi 

vE 
tE e 
-pE 



zE 
wE\ 
xE 

-uEy 

-zE 
wEi 
xE 
uEi 



pE 
yEi 
vE 
-tE x 
-pE 
yEx 
vE a 



xE 

uEx 

-zE 

-wEx 

xE 

-uEx 

zE 

-wEx 



yEx 

vE Q 
tEx 
-pE 
yEx 
vE 
-tEx 
pE 



-zE 
wEx 
xE 
uEx 
zE 
wEx 
xE 
-uEx 



-tEx \ 

- P E \ 

yEx 

vE 

tE t 

pE 

yEx 

vE 



Then W is a paraunitary matrix in the variables {x, y, z, t, u, v, w,p}. 

Interchanging blocks of rows and/or columns in W will also give a paraunitary matrix which is 
equivalent to W. By varying the signs other constructions of paraunitary matrices are obtained and 
these are not generally equivalent to one another. 
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They may seem to be non-separable but by interchanging blocks of rows and blocks of columns it can 
be shown that they are separable. However they are interesting in themselves with interesting properties 
and can also be used to construct Hadamard matrices of regular types. 

Examples may also be interpreted as coming from the structure of group rings of various groups. 
Then giving the values ±1 to the variables can result in regular Hadamard matrices or giving the values 
e lB to the variables can result in Hadamard complex matrices. Hadamard matrices with entries which 
are roots of unity may also be obtained from these constructions. 

The following is an example of this type. 



w(x,y,z,t) = - 



v -y 
-y y 



y -y 
-y y 



-t t 



-y y z 



y -y 
-y y 



y -y 
-y y 



-y y 
y -y 



t -t 
-t t 



i -v y 

\ v -y 



-y y 

y -y 



-t t 



-t t 



-t t 



y -y 
-y y 



-t 
t 



y -y 

-y y 



where PQ* = = QP* and has the structure of the group ring 



This matrix has the form 
of C 8 x C 2 . 

It is clear that x, y, z, t may each be replaced by a monomial times a complex number of modulus 1 
in W and a paraunitary matrix is obtained. 

Values of ±1 may be given to the variables in W and with the fraction omitted this gives a regular 
real Hadamard matrix. Values of modulus 1 may be given to to the variables in W and with the fraction 
omitted a regular complex Hadamard matrix is obtained. Other group ring structures can also arise in 
this manner. 

Walsh-Hadamard matrices have the structure of the group ring of Cg- These examples may also be 
extended in a similar way. So for example constructions with the structure of the group ring C% x C2 x G 



may be made if paraunitary matrices with the structure of G can be formed such as when G 
when G = (C 2 x C 2 ) n and others. 



C? or 



7 Non-separable constructions 

Several methods have now been developed for constructing multidimensional paraunitary methods from 
complete orthogonal sets of matrices. The matrices produced are useful in many ways but are found to 
be separable although not appearing so initially. The methods use just one complete orthogonal set of 
idempotents in the various constructions. 

Thus we are led to consider different complete orthogonal sets of variables and to 'tangle' them up 
in order to construct non-separable paraunitary matrices. 

7.1 A general construction 

Proposition 7.1 Let A, B be paraunitary matrices of the same size but not necessarily with the same 
variables over a field in which 2 has a square root. Then W = ^^j ^B^j an ^ ® = \73 
are paraunitary matrices in the union of the variables in A,B. 
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Proof: Suppose A, B are n X n matrices. Then 

ww* = l -( A B )( A * A * 

2\A —B J \B* —B* 

1 (AA* + BB* AA* - BB" 

2 \AA* - BB* AA* + BB" 

_ 1 (In In In Ir. 
2 V n In In ~\~ Ir. 

The case for Q can be considered similarly or it follows from Lemma 13-11 since Q is the transpose of 
W. □ 

Paraunitary matrices constructed by methods of previous sections may be used as input to Proposition 
IT. II to construct paraunitary matrices. Matrices constructed using the Proposition can then also be used 
as input. 

The methods are fairly general and it is easy to produce examples for input using various complete 
orthogonal sets of idempotents. The result holds in general over any field which contains the square root 
of 2. 

If A = B then W in Proposition 17.11 is the tensor product A® J where J — ( } 1 1 ) . If A and B are 
formed using the same complete symmetric orthogonal set of idempotent as in [5] or H2 then W can be 
shown to be separable. 

I vA aB 

It would appear initially that Proposition 17.11 could / should be generalised to W — -j= [ A ^ 

where p, q are monomials or in the case of C where p, q are monomials times a complex number of 

f A B \ (pi 



modulus 1. However then W is separable as a product W — ■ ^ 

If W = ( „ ^p] where X,Y,Z,T are matrices of the same size then X,Y,Z,T are referred to as 



, Z T , 

the blocks of W and [X Y) and {Z T) are the row blocks of W . Similarly column blocks of W are 
defined. 

Suppose A, B are matrices of the same size. Then a tangle of {^4, B} is one of 
A B 



l - W -^\A -B 

2. A matrix obtained from 1. by interchanging rows of blocks and/or columns of blocks. 

3. The transpose of any matrix obtained in 1. or 2. 

A tangle of {A, B} is not the same as, and is not necessarily equivalent to, a tangle of {B,A}. 
Note that interchanging any rows and/or columns of a paraunitary matrix results in an (equivalent) 
paraunitary matrix. Thus in particular interchanging rows and/or columns of blocks also results in 
equivalent paraunitary matrices; thus item 2. gives equivalent paraunitary matrices to item 1. The 
negative of a paraunitary matrix is a paraunitary matrix as is the * of a paraunitary matrix. 

For example 

75 (1 - B ) ' 71 (1 "/) ' 75 (s - B ) are tangles of {A > B} 
and 

75 (b -a) ' 75 [b ~a) ' 75 [a -a) are tangles of {B > A} ' 
Proposition 17. II may be generalised as follows. 



Proposition 7.2 Let A, B be paraunitary matrices of the same size but not necessarily with the same 
variables. Then a tangle of {A, B} or {B, A} is a paraunitary matrix. 

Use the expression 'A is tangled with B^ to mean that a tangle of {A, B} or {B, A} is formed. 
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7.2 Examples 

1. (a) Construct A — (x) and B = (y). 

(b) Construct W = ( ) = ( % - y ) . Then is a paraunitary matrix. 

(c) Similarly construct Q = ( \ f t ) . 

(d) Tangle W and Q to produce for example the paraunitary matrix T = 



x y z t 
x —y z —t 
x y —z —t 
jc—y—z y 



(e) The process can be continued: Matrices produced from (d), with different variables, can be 
input to form further paraunitary matrices. 

(a) Construct as is 14.81 the following complete symmetric sets of idempotents in 3 x 3 matrices 
over F7: 

^=(!H).« = aiI).^(If!)>.w»=(!Ii).«-(IlI)^=(§8!)>- 

(b) Form A = xP + yP x + zP 2 ,B = tQ + rQi + sQ 2 . 

(c) Tangle A, B to form for example the following paraunitary matrix over F7: ( J 



A A ' 
B B , 



(a) Construct, in C2X2, the following complete symmetric (different) sets of orthogonal idempo- 
tents {E ,Ei} and {Qo,Qi} where: 

£ = |(ll):- E l = H-l~l 1 )' Qo = 5 ( I 1 ) ' Ql = 5 ( -2 ~i ) ■ 

(b) Construct A — xE + yEi,B — zQ + tQ\. 

(c) Construct W = -j= ^ ^B^j ' ^ en W is a. paraunitary matrix of size 4x4 with variables 
{x,y,z,t}. 

(a) Construct different {Eg, E{\ and {Qd Qi} complete symmetric orthogonal sets of idempotents 

in C nX n- 

f xE Q yEi uQ vQi ^ 
rEi pE zQi tQ 
xE yEi -uQq -vQ\ 
\rEt pE -zQi -tQoJ 
It is essential that {Eq, Ei} and {Qq, Qi} are different complete orthogonal sets of idempotents 
in order for W to be non-separable although the construction does not depend on this. 

(c) Clearly also the roles of {Eq, Ei} and {Qo,Qi} can be interchanged in W and a paraunitary 
matrix is still obtained. Changing the ± signs in such a way that the block inner product 
of any two rows of blocks is will give a different inequivalent paraunitary matrix. Thus for 
example W could be replaced by the following: 

/ xE yEi uQ vQi \ 

rEi pE -zQi —tQo 

-xE Q —yEi uQ vQi 

\ rE x pE Q zQi tQ J 



(b) Construct W = -j- 



Then W is a paraunitary matrix. 



(d) Construct W = 4= 



Then this W is also a paraunitary 



matrix. 



(a) See section|4]for methods for constructing complete orthogonal symmetric sets of idempotents. 
Examples of such sets in C2X2 are {Eq, Ef}, {Qq, Q\} where these are given as in Example 2 
above. 

(b) Another complete symmetric orthogonal set of idempotents in C2X2 is the following: 
{ p o = \ (] 7) ' p i = 2 (-M)}' 

(c) In example 2. {Eq,Ei} is 'tangled' with {Qo,Qi}- {Po,Pi} may similarly be combined 
('tangled') with either {Eq,Ei} or {Qo,Qi} to construct paraunitary matrices. 

(d) Using {Po, Pi} with {Eo, Ef} produces paraunitary matrices of the form ^P where the entries 
of P are ±l,±i with i = y/—l. By specialising the variables, complex Hadamard matrices 
may be obtained. 



16 



x x y —y u 



/ x x 

X X 



-IU V 



7 V 



(c) For example the following is a paraunitary matrix: 



-it 



r —r p p z iz t 
— r r p p —iz z it t 
x x y —y —u iu —v —iv 
x x —y y —iu—u iv —v 
r p p —z—iz—v iv , 
— iz z ' 



r r p p iz 



(f) By giving values of modulus 1 to the variables, complex Hadamard matrices are obtained. 

For example letting all the variables have the value +1 gives the following complex Hadamard 

matrix: 

/ 1 1 1 -1 1 — i 1 i \ 

1 1-11 i 1 -i 1 \ 

1-11 1 1 i 1 -i 

-11 1 1 -i 1 t 1 
1 1 1 -1 -1 i -1 -i 

1 1 -1 1 -i -1 i -1 



\ 1-11 1 -1 -i -1 i i 
\ -1 1 1 1 i -1 -i 1 / 

6. (a) Construct Pq = 5(212^, Pi = | 



12-2 
2 4-4 
-2 -4 4 



,P> 



4 -4 -2 
-4 4 2 
-2 2 1 



(b) Construct the complete symmetric orthogonal set of idempotents obtained from the group ring 
CC3 of the cyclic group of order 3: Qo = circ(l, 1, 1), Qi = circ(l, u, oj 2 ), Q 2 = circ(l, uj 2 , v) 
where a; is a primitive cube root of 1; 

^tFo yPi zP 2 \ (ay 

dQ 2 

K gQi 



(c) Construct A 



P P 2 qP rPi 
,sP 1 tP 2 vP 0/ 



and B 



bQi 

eQo 
hQ 2 



(d) Construct W = ^ ( A . 



Then W is a paraunitary matrix. It has size 18 x 18 and 18 variables; variables can be equated. 



W in the above could for example be replaced by W 



1 

V2 



where p, q are variables and 



' ipA qB 
ipA —qB 

i = \J — 1 but as pointed out this is separable and may be constructed as a product. 

A complex Hadamard matrix is a matrix H of size n x n with entries of modulus 1 and satisfying 
HH* = nl n . Complex Hadamard matrices arise in the study of operator algebras and the theory of 
quantum computation. 

By giving values which are k th of unity to the variables, with k divisible by 4, in the above example 
4, special types of complex Hadamard matrices which are called Butson-type are obtained. A Butson 
type Hadamard H (q, n) matrix is a complex Hadamard matrix of size n x n all of whose entries are q th 
roots of unity. 



Here is another example which uses group rings:. 

1. Construct the complete symmetric set of orthogonal idempotents {Pi\i — 0,1,..., 5} from the 
group ring CCq of the cyclic group Cq of order 6. This gives P,; = circ(l, a/, uj 21 , w 3j , lu 4 ' 1 , lu 5 ' 1 ) 
where w is a primitive 6th root of unity. 

2. Define Q = Po,Qi = P\ + P5, Q2 = P2 + Pi- Note that Q , Qi,Q 2 are real. 

3. Let 



niiir 
111111 
111111 
111111 
111111 
.111111, 



,E 2 = l 
6 



1-1-1-11 1 
-11 1 1-1-1 
-11 1 1-1-1 
-11 1 1-1-1 
1-1-1-11 1 
1-1-1-11 1 



,Ea = — 



2 -1 -1 - 
2-1-10 
0-12-10 
0-1-12 
-10 2-1 
.-1 0-12 



be the complete symmetric set of orthogonal idempotents obtained from the group ring of Ss(= Dq) 
as in section 14.71 



' xEq yEi zE 2 • 
4. Define A = I pE 2 qE rE x | and B 

sEi tE 2 vE , 



f aQa bQi cQ 2 
dQ 2 eQ fQi 
K gQx hQ 2 kQ 0/ 



17 



(a b 

5. Construct W — -j- ( ^ ^ 

Then W is a paraunitary matrix which is real. It is a 36 x 36 matrix with 18 variables. 



3 An Algorithm 

1. Construct different sets {Pq, P2, ■ ■ ■ , Pk} and {Qo, Qi, . . . , Qk} of complete orthogonal symmetric 
of idempotents by the methods of section I4TT1 or section I4T31 for otherwise) in F nxn . (F is usually 
C but other fields can also be used.) 

2. Construct a paraunitary matrix from {Pq,P2, ■ • • , Pk} by either the methods of section [5] or the 
methods of section [HI Call this matrix A. 

3. Construct a paraunitary matrix from {Qo, Qi, . . . , Qk} by either the methods of section [5] or the 
methods of section [6j Call this matrix B. The variables in A, B can be different. 

4. Construct a tangle of {A, B} or {B,A}. 



An example: 



1. Construct {P , Pi, P2, P3} and {Qo, Qi, Q2, Q3} in C4 X 4 where the Pi and Qj are obtained by 
construction methods of 14.11 14. 



2. Form 



/Pa Pi Pi P 3 \ 

P3 Pq Pi P2 

P% Pz Pq Pi 

\Pt P 2 P 3 P J 



(Here the structure of C4 is used and a circulant structure is obtained.) 



3. Form 



x 3 iP 3 \ 
x 3 iPi 
x 33 P\ 

XuPqJ 

Let the matrix in 3. be denoted by A. 



(x m P XnPi X21P2 

x 2P 3 X12P0 X22P1 

X 03 P2 x 13 P 3 x 23 Pq 

\x04P1 x 14 P 2 X24-P3 



5. Form 



6. Form 



/Qo 
Qi 
Q 2 

\Q 3 



Qi Q 2 

Qo Q3 

Q3 Qo 

Q 2 Qi 



Q 3 \ 
Qi 
Qi 
QoJ 



(Here the structure of C2 X C2 is used.) 



/z/oiQo VuQi 

2/02Q1 yi2Qo 

2/03Q2 2/13Q3 

V2/04Q3 2/14Q2 



X21Q2 
ynQ 3 
2/23Q0 

2/24 Ql 



a^3iQ3\ 
Z/32Q2 

2/33 Ql 
2/34 Qo / 



7. Let the matrix in 6. be denoted by B. 

8. Form W = 



B 

-B 



4 Further algorithm 

1. Input paraunitary matrices A, B of the same size but not necessarily with the same variables. 
These may be formed from method of section 17.31 or from this algorithm. 

2. Form a tangled product of {A, B} or {B, A}. 
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7.5 Further constructions 

The non-separable paraunitary matrices and separable paraunitary matrices can be combined when 
appropriate as products or as tensor products to construct further paraunitary matrices. These may 
then be input to algorithm of section 17.41 

8 Pseudo-paraunitary 

Let P be a paraunitary n x n matrix with variables z over a field F. Then the rows {v%, v 2 , ■ ■ ■ , v n } of P 
satisfy v^* — 1 and ViV* = for i ^ j. Note that v* means transpose conjugate over C, and transpose 
over other fields, with the understanding that z* = z ,{z }* = z for a variable z. 

Let Pi — v*Vi which are n x n matrices of rank 1 and involve the variables {z, z -1 }. Then 
{Pi, P 2 , . . . , P n } is a complete orthogonal symmetric set of idempotents in the polynomial ring F nxn [z, z _1 ] 
Hence by the methods of Sections01ini[S]and[7]paraunitary-type matrices may be formed; for the method 
of Section[7]two such sets must be constructed. For example W = w\P\ + w 2 P 2 + . . . + w n P n in variables 
w = (wi, w 2 , ■ ■ ■ , w n ) satisfies WW* = 1. Now W is a matrix in the variables (z, z _1 , w) but cannot be 
termed paraunitary. Call such a matrix a pseudo-paraunitary matrix. Having constructed W its rows 
may then be used to construct further pseudo-paraunitary matrices and so on. 

Thus say W(z, z _1 ) G F nxn [z, z" 1 ] is a pseudo-paraunitary matrix if WW* — 1. Pseudo-paraunitary 
matrices may be constructed from paraunitary matrices and from pseudo-paraunitary matrices. 

Here's an example. 

1. ¥ovmE = \{\\),E l = \{} 1 - 1 1 ). 

2. Form P = xE +yE 1 = 1^ 

3. Let vi = i(x + y,x-y),v 2 = \{x-y,x + y). 

4. Form Pi = v*v\ = 4 , y u , Po = v*v 2 = 4 _i _i _i -1 • 

1 11 4 ^ x x y~y 1 x 2-x y-y • x / ' 1 4^2; 1 y-y x x 2+x x y+y l x J 

5. pp = p, PiP 2 = and Pi + P 2 = 1. Form PF = zP x + tP 2 . 

6. W = W{x 1 y,x- 1 ,y- 1 ,z,t) 1 and WW* = 1. 

7. The rows of VF can be used to form further pseudo-paraunitary matrices. 

Consider Pi,p2 as in this example above. Define Q\ — xyP\,Q 2 = xyP 2 . Define Q = zQ\ + tQ 2 - 
Then Q 6 P2x2[^, y, x, t], and is a polynomial with QQ* — x 2 y 2 I 2 - 

Say W £ P nX n[z] is a pseudo-paraunitary matrix over P nx „[z] if WW* = pl n where p is a monomial 
in z. A pseudo-paraunitary matrix over P„ x „[z,z _1 ] may be used to construct a pseudo-paraunitary 
matrix over P„ xn [z] and vice versa. 

Pseudo-paraunitary matrices in general may be constructed from paraunitary matrices and from 
pseudo-paraunitary matrices. 

9 Determinants and rank 

Here we consider properties of complete sets of idempotent matrices and ranks of the idempotents. 

Lemma 9.1 Suppose {Pi, E 2 , . . . , E s } is a set of orthogonal idempotent matrices. Then rank(Pi + E 2 + 
... + E s ) = tr(E x + E 2 + . . . + E s ) = trEi + trE 2 + ... + trE s = rank E x + rank E 2 + . . . + rank E s . 

Proof: It is known that rank A = tvA for an idempotent matrix, see for example [2], and so rankP^ = 
trEi for each i. If {E,F,G} is a set an orthogonal idempotent matrices so is {E + F.G}. From 
this it follows that rank(Pi + E 2 + . . . + E s ) = tr(Pi + E 2 + . . . E s ) = trPi + trP 2 + • • ■ + trP s = 
rank Pi + rank E 2 + . . . rank E s . □ 
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Corollary 9.1 rank(£' il + E i2 + . . . + E ik ) = rank £7^ + rank£"j 2 + . . . + rank_E i(c for ij e {1,2, ... , s}, 
and ij ^ it for j ^ I. 

Let {ei, e 2 , . . . , e k } be a complete orthogonal set of idempotents in a vector space over F. 

Theorem 9.1 Let w — a.\e\ + a 2 e 2 + . . . + a k e k with on e F. Then w is invertible if and only if each 
a; + and in this case w^ 1 — —e\ + — e 2 + . . . + — e k . 

Proof: Suppose each a l ^ 0. Then w(^e + ^-ei + . . . + -±-e k ) = e§ + ef + . . . + e\ = e + ei + . . . + e fc = 1. 

Suppose w is invertible and that some Qj = 0. Then me^ = and so to is a (non-zero) zero-divisor 
and is not invertible. □ 

We now specialise the Ci to be n x n matrices and in this case use capital letters and let e; t = Ei. 

Let A = a\E\ + a 2 E 2 + ■ ■ ■ + a k E k . Then A is invertible if and only if each ai ^ and in this case 
A- 1 = —Ei + ±E 2 + ... + ±E k . 

ai 1 a 2 z a k K 

Theorem 9.2 Suppose Ei,E 2 , ...,E k is a complete symmetric orthogonal set of idempotents in F nxn . 
Let A = a 1 E l +a 2 E 2 +. . .+a k E k with a t e F. Then the determinant of A is \A\ = El a r 2 ank E2 . . . a k ankEk . 

Proof: Now AEi — aiEf — aiEi. Thus each column of Ei is an eigenvector of A corresponding to 
the eigenvalue a^. Thus there are at exist ranki?, linearly independent eigenvectors corresponding to 
the eigenvalue ai. Since rank£?i + rank£?2 + ... + r&nkE k = n, there are exactly ranki?i linearly 
independent eigenvectors corresponding to the eigenvalue ai. Let = rank^. Let these Ti linearly 
independent eigenvectors corresponding to dj be denoted by Vi,2, • • • v», n . Do this for each i. 
Any column of Ei is perpendicular to any column of Ej for i ^ j as EiE* = 0. 

Suppose nOW Y^j=l <*l,jVl >rj + YJjLl a 2j W 2,r 3 + • ■ • + X^li a k,j = 0- 

Multiply through by E s for 1 < s < k. This gives Y^j=i a k,jV k j = from which it follows that 
a k ,j = for j = 1,2,. ..r k . 

Thus the set of vectors S = {«i,i,ui,2, • • • i>i,n , «2,i, v 2 ,2, • • • , «2,r 2 ,Ufe,i,«fc,2, • • - ,v k ^ k } is lin- 

early independent and form a basis for F n - remember that rank(£?i + E 2 + . . . + E k ) = n. Hence A 
can be diagonalised by the matrix of these vectors and thus there is a non-singular matrix P such that 
P~ 1 AP = D where D is a diagonal matrix consisting of the at repeated times for each i — 1, 2, ... ft. 

Hence \A\ = \D\ =a r 1 1 a r 2 2 ...a r k k . □ 
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